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Abstract 

We present a Bethe approximation to study lattice models of linear poly- 
mers. The approach is variational in nature and based on the cluster variation 
method (CVM). We focus on a model with (i) a nearest neighbor attractive 
energy e v between pair of non-bonded monomers, (u) a bending energy for 
each pair of successive chain segments which are not collinear. We determine 
the phase diagram of the system as a function of the reduced temperature 
t = — and of the parameter x = — . We find two different qualitative beh&v- 
iors, on varying t. For small values of x the system undergoes a 6 collapse 
from an extended coil to a compact globule; subsequently, on decreasing fur- 
ther i, there is a first order transition to an anisotropic phase, characterized 
by global orientational order. For sufficiently large values of x, instead, there 
is directly a first order transition from the coil to the orientational ordered 
phase. Our results are in good agreement with previous Monte Carlo sim- 
ulations and contradict in some aspects mean-field theory. In the limit of 
Hamiltonian walks, our approximation recovers results of the Flory-Huggins 
theory for polymer melting. 

PACS numbers: 05.70.Fh, 36.20.Ey, 64.60. Cn 

The configurational statistics of a long, linear polymer in solution has often been modeled 
by a self-avoiding walk (SAW) on a lattice |],@. The self-avoiding constraint takes into 
account excluded-volume effects. Attractive Van der Waals interactions between monomers 
are also generally considered. They are included by assigning a negative energy — e v to each 
pair of nearest-neighbor monomers on the lattice, provided they are not consecutive along 
the chain. These interactions become relevant at low enough temperature, causing a collapse 
transition of the polymer ||. The transition point is called 9 point and it separates a high 
temperature expanded structure from a low temperature compact globule. 

Polymers with a local stiffness have also been considered, by including a bending energy 
th which favors straight segments of the chain &-H. The semi-flexible model has attracted 
much interest in the low temperature phase and in particular in the limit of Hamiltonian 
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walks, where the path is forced to visit all sites of the lattice. In this case, it is believed 
to describe the melting of polymers chains 0] (see also 0). The system undergoes a phase 
transition between a disordered (liquid) phase and an ordered (solid) phase, the latter being 
characterized by anisotropic orientational order. More recently, the semi-flexible model 
has attracted renewed interest because of the possible connection with the protein folding 
problem [§,0. In this spirit, each link of the walk represents an a-helical turn (ca. 3 amino 
acids) and the curvature term mimics the tendency to form secondary structures. The 
attractive energy between monomers models the hydrophobic effect, which is supposed to 
be the main driving force for the folding transition ]TU| . 

In the present work we consider a Bethe approximation for lattice homopolymers. Our 
aim is twofold: first we introduce a new method to deal with SAW problems, second we 
present results concerning the phase diagram of an isolated semi-flexible chain. In the 
approach to the problem we have followed the cluster variation method (CVM) pTl|— jl3|| . 
This is a closed form approximation, which is known to give excellent results for the phase 



diagram of spin systems |14| . The approximation scheme is determined by the largest clusters 
of sites which are treated exactly. The CVM allows us to write an approximate expression for 
the free energy of the system, as a function of the probability of occurrence of all possible 
configurations of the basic cluster. This free energy has then to be minimized, subject 
to consistency conditions on the distribution variables. The pair approximation considers 
nearest-neighbor pair of lattice sites and it corresponds to the Bethe approximation. 

We represent a polymer as a SAW on a <i-dimensional hyper-cubic lattice with V = 
L d sites. Pairs of non-consecutive vertices along the chain interact through an attractive 
nearest-neighbor potential e v . Stiffness of the chain is incorporated by attributing an energy 
penalty to each turn (corner) of the walk. Let T be the absolute temperature and (3 = -^Kp 

the inverse temperature. In the following we will adopt the notation uj = f3e v , t = — and 
x _ it partition function of the system is 

Z N = V" e uj(N con (S)-xN cor (S)) ^ 
{SAW} 

where denotes the ensemble of all iV-step SAW; N con (S) and N cor (S) are respec- 

tively the number of contacts and corners in walk S. Introducing a monomer fugacity z, the 
grand canonical partition function reads therefore 

oo 

Z = £ z N Z N (2) 

N=l 

where the sum is over all possible lengths N of the walk. 

The lattice model ([!]) (or equivalently (fj)) has been the object of recent investigations 
0-^. A mean-field analysis predicts a ^-collapse transition at a temperature t g inde- 
pendent of x. Another transition should occur at lower temperature. It is a discontinuous 
melting transition from a disordered globule to an ordered "crystalline" phase. The melting 
temperature t m increases with x, although in mean-field approximation t m < t e for any value 
of x. This picture has been partly contradicted by heuristic arguments HH and by Monte 
Carlo simulations f7],§]. Indeed t e appears to slightly increase with x. Most importantly, t m 
seems to grow with stiffness and it does not reach a finite asymptotic value. This implies 



2 



that the line of the melting transition hits the line of the 0-collapse. For sufficiently high 
values of x, there is therefore a direct first order transition from the open coil to the ordered 
phase. The triple point is found to be approximately at x — 13 M. 

The formulation of the CVM given by An [13|] is a particularly convenient starting point 



for our analysis. We first illustrate the method by treating in some detail the simpler case 
of zero stiffness (x — 0). Let Sj and pi be the distribution variables assigned respectively 
to each site and pair configurations. One should in principle distinguish among all the 
possible configurations which are not related by symmetry operations. In d — 2 there 
are, for instance, 3 single site and 11 pair independent configurations. In fact we have 
verified that in the particular case x = they can be grouped into a smaller number of 
non-equivalent classes. These are determined only by the following conditions: (a) the 
site is visited by the path, (b) the nearest-neighbor pair is joined directly by the path. 
The independent configurations are reported schematically in fig. |I|, together with their 
multiplicity of occurrence. Following the notation of fig. Ill the free energy of the system 



in the pair approximation, can be written [13 



(3F q(q-l) , g(g~l) 2 (g-2) 2 

— = In zsi Lup 2 + 




In s^j + | m p(S)Vi ^Pij (3) 



where q = 2d is the coordination number of the lattice and m s (i) and m p (i) stand respectively 
for the multiplicity of site and pair configurations. Normalization of the distributions and 
consistency conditions on the probability variables require respectively 

S2 = 1 _^ Sl (4 ) 



and 

Pi 



P2 " ( q -l)(q-2) - (q-l)(q-2f 3 ^ 

q(q-l) (g-lXg-2) 
P4 = 1 - —2—*i 2 P3 

This leaves us with only two variational parameters, e.g. s\ and p%. Substituting (Q) and 
(^) into ([3]) and minimizing with respect to s% and p^ we obtain the stable phase at a given 
z and io. We report the complete phase diagram for d = 3 in fig. |2]. The polymer is a critical 
system along the transition line z c (u). This line separates a chain with zero density (si = 
for z < z c (u)) from a chain with finite density (s\ ^ for z > z c (u)). The continuous 
line represents a second order transition and the average number of monomers diverges with 
a power law as z tends to z c . The broken line is instead a first order transition and the 
density of monomers makes a finite jump at z c . The cross denotes the tricritical point and it 
corresponds to the 9 point |Tj| . In the case of pure SAW (uj = 0) the connectivity constant is 



_1 = 2c? — 1. This result would have been expected by studying SAW on a Bethe lattice 
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and it should be compared, for instance, with /i rs 2.64 |Tg] and /x ~ 4.68 flT7 
from exact enumerations respectively in d = 2 and d = 3. In our framework z, 



obtained 
does not 



depend onw, as long as u < to g . This is certainly an artifact of the approximation. It can be 
ascribed to the fact that, after minimization in p 3 , p 2 ~ in the limit of si going to zero. 
As a consequence, there is no term in ([3]) proportional to uj which contributes in locating 
the minimum of the free energy around s\ — 0. Nonetheless the estimates we obtain for the 
9 point uj { b b) are a better approximation to the available numerical values, with respect to 
mean field theory (u ( e MF) = ± §): in d = 2, u e « 0.665 @, « 0.4055 and uj g MF) = 0.25; 
in d = 3, u B » 0.275 0, c<4 s) « 0.2231 and cj< Mi ^ = 0.1667. 

In the general case x ^ equations (0), (|j) and (^) must be generalized to include the 
curvature energy and the possibility of an anisotropic phase. In this case it is not possible 
to group configurations as in fig. [I] and one has to face a complex minimization problem. An 
efficient way of doing it numerically is by mean of the natural iteration method |12| . The 
resulting phase diagram is reported in fig.[| as a function of x and t for d = 3. The fugacity 
z is fixed to its critical value z c (x,t). This condition assures we are studying a polymer in 
the limit of infinite chain length (N — > oo). We find three different phases: an open coil, a 
compact globule and an ordered crystal. In our approximation, the latter is just the ground 
state of the polymer, having all links perfectly aligned. This is known to be not completely 
correct, as it has been shown rigorously that, for instance, in the case of Hamiltonian walks 
the entropy strictly vanishes only in the limit T —>■ [^(],^l[] . The ^-collapse line between the 



coil and the globule appears to be independent of x. On the other hand the discontinuous 
melting transition tends to infinity with x. Beyond the triple point at x ~ 8.8, there is 
directly a first order transition from the coil to the solid. 

The limit T — > (or z — > oo) corresponds to Hamiltonian walks. In this limit, walks are 
space filling and configurations with vacancies do not contribute to the partition function (0). 
Also, the attractive nearest-neighbor potential e v plays no role as there are precisely (d — 1) 
contacts per monomer. The resulting model is the so called Flory model of polymer melting 
II. In this case, we obtain a minimum for the free energy in the disordered, compact phase 
which coincides with the free energy estimated by using the Flory-Huggins approximation 
22|. The latter was originally derived from combinatorial arguments ||23[| . Neglecting the 



constant contribution proportional to e v , the analytical expression reads 



(3F F 



V 



In 



1 + (q — 2) exp 



-LUX] 



(6) 



At low temperature F FH competes for stability with the local minimum associated to the 



ordered phase, which has strictly F D 
takes therefore place at 



in our approach. A first order phase transition 



x 



In 



g-g 



'-i 



(7) 



In particular for d = 3 (q = 6) we have t m /x = fin (J§-Y) ~ 0.86, which corresponds to the 
slope of the globule-solid transition line of fig. in the limit x — > 0. This value is slightly 
larger than the analogous mean-field estimate, t m /x ~ 0.58 [25 . 
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In the case x = each HW is equally weighted and the total number of paths is believed 
to scale as M HW — //^. From @ we therefore obtain 



A mean-field approach |[24|| , in very good agreement with numerical estimates, yields 



/4 MF) = \ (9) 

It is interesting to note that expression (|8]) correctly predicts /j, h = 1 for d = 1, contrary 
to @. In order to evaluate corrections to mean-field theory in powers of 1/q, we have 
formulated the Hamiltonian walk problem through spin variables. This is in close analogy 
to the De Gennes theorem for SAW [p6fl . By using a suitable high-temperature expansion 
we have then been able to calculate the coefficients of the series up to the third order, 



extending of one order previous results |28|. We find 



"» = ; (1 + 6?-3? + -> < 10) 



In the limit q — ► oo expression (H) specializes to 



The Bethe approximation therefore correctly recovers results of the expansion (jT0|) up to 
order O(^). 

In summary, we have presented a Bethe approximation for lattice models of linear ho- 
mopolymers. The method constitutes a substantial improvement with respect to mean-field 
theory. Indeed it produces a phase diagram for a semi-flexible polymer chain which is in 
good agreement with Monte Carlo simulations. In particular, we find a triple point where 
the 9 collapse line and the melting transition line meets. In the limit of Hamiltonian walks it 
recovers results of the Flory-Huggins theory for polymer melting, whose variational nature 
appears in a transparent shape within our framework. It has the advantage of not requiring 
any spin or field theoretical representation, rather it relies directly on the configurations of 
the system. This last consideration suggests the scheme is more general and suitable to be 
applied to other geometrical problems, as, for instance, branched polymers p9f and self- 



avoiding surfaces [^] (see and references therein). It is plausible that the accuracy of the 
method can be systematically refined according to the cluster variation method, in analogy 
with spin systems. We expect some of the inaccurate features of the pair approximation 
could be removed by considering larger basic clusters. 

a Electronic address: lise@sissa.it 

b Electronic address: maritan@sissa.it 

c Electronic address: alex@athena.polito.it 
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FIG. 1. Schematic representation of independent (a) site and (b) pair configurations, in the 
case of zero stiffness (x = 0). The continuous line represents the path visiting a site; q = 2d is the 
coordination number of the lattice. 
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FIG. 2. Phase diagram of the system as a function of lo and z, in the case of zero stiffness 
(x = 0). The average length of the polymer is finite (infinite) in region I (II). The continuous 
(dashed) line is a second (first) order transition. The cross marks the tricritical point (to g ~ 0.2231 
and z g = 0.2). 




FIG. 3. Phase diagram of the system as a function of x and t. The continuous line denotes 
the 9 transition from the coil (C) to the globule (G). The dashed line represents the first order 
transition to the solid (S). The triple point is at x ~ 8.8 and t « 4.5. See fig. 3 in Q and fig. 8 in 
H for comparison with Monte Carlo simulations. 
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